Gaussianizing the non-Gaussian lensing convergence field II: the applicability to noisy 
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In paper I (Yu et a/. we showed through N-body simulation that a local monotonic Gaus- 
sian transformation can significantly reduce non-Gaussianity in a noise-free lensing convergence 
field. This makes the Gaussianization a promising theoretical tool to understand high-order lensing 
statistics. Here we present a study of its applicability in lensing data analysis, in particular when 
shape measurement noise is presented in lensing convergence maps, (i) We find that shape measure- 
ment noise significantly degrades the Gaussianization performance and the degradation increases 
for shallower surveys, (ii) The Wiener filter is efficient in reducing the impact of shape measure- 
ment noise. The Gaussianization of the Wiener-filtered lensing maps is able to suppress skewness, 
kurtosis, and the 5th- and 6th-order cumulants by a factor of fO or more. It also works efficiently 
to reduce the bispectrum to zero. 

PACS numbers: 98.80.-k; 98.65.Dx; 98.62.Ve; 98.62.Sb 



o 

6 



> 
o 



I. INTRODUCTION 

Weak gravitational lensing — a powerful probe of the 
dark universe — is facing many challenges. One of them 
is the nonlinear evolution of the large scale structure. 
The non-Gaussianity induced by the nonlinear evolu- 
tion pushes the cosmological information into high-order 
statistics, leading to the loss of constraining power from 
the 2-point statistics (the power spectra and equivalently 
the correlation functions) alone. Recently several works 
have focused on nonlinear local transformation of the 2D 
lensing and 3D matter and galaxy fields to reduce the 
non-Gaussianity [H-Ql- This procedure is shown to be 
promising. The nonlinear transformation changes the in- 
formation distribution in the hierarchical statistics. It 
makes the cosmic fields close to Gaussian and hence sig- 
nificantly improves the information we can extract from 
low-order statistics. 

Most of these works adopt the logarithmic transfor- 
mation. Instead, we choose a Gaussian transformation, 
which is defined as a local monotonic transformation to 
Gaussianize the one- point PDF Throughout the pa- 
per we refer to the application of this transformation as 
the Gaussianization. This approach is motivated by the 
Gaussian copula hypothesis [8] . If this hypothesis is valid, 
the Gaussianization will make all the statistics Gaussian, 
in addition to the one-point PDF. 

Through N-body simulations, we quantified the per- 
formance of the Gaussianization on a 2D lensing field 
by various measures of non-Gaussianity, such as skew- 
ness, kurtosis, the 5th- and 6th-order cumulants as func- 
tions of smoothing radius, and the bispectrum. We found 
that the Gaussianization works surprisingly well. It sup- 
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presses the above non-Gaussianity measures by orders of 
magnitude and effectively reduces them to zero. This im- 
plies that in many exercises we can treat the new field 
as Gaussian. The Gaussianization procedure then trans- 
ports cosmological information into the power spectrum 
of the new one. For this reason, theoretical modeling of 
the weak lensing statistics can be significantly simplified. 

However, this does not necessarily mean that the Gaus- 
sianization procedure works for real data, in which there 
are various measurement noises. A prominent one is 
the shape measurement noise from the random galaxy 
shapes. Its impact on the Gaussianization can be figured 
out straightforwardly in two limits: (i) When the shape 
noise is negligible with respect to the lensing signal, the 
Gaussianization works rather well, as shown in paper I. 
(ii) When the shape noise overwhelms the lensing signal, 
the Gaussianization collapses. Let us look at the case 
of a sufficiently large pixel with ^ 1 galaxies. Due 
to the central limit theorem, the shape noise distribu- 
tion approaches a Gaussian, as does the one-point PDF. 
Hence the Gaussian transformation becomes trivial and 
the Gaussianization fails. 

Hence, whether or not the Gaussianization is applica- 
ble to lensing surveys is a nontrivial question. To inves- 
tigate the applicability, we add shape measurement noise 
to the simulated lensing convergence maps. Following 
paper I, we then test the performance of the Gaussian- 
ization through several measures of non-Gaussianity, in- 
cluding the cumulants up to 6th-order of the smoothed 
field, and the bispectrum. 

We find that the performance of the Gaussianization is 
strongly affected by shape measurement noise. However, 
if we Wiener filter the lensing maps first and then perform 
the Gaussianization, it can still work rather well for DES- 
like or LSST-like surveys. 

Our paper is organized as follows: In Sec|lT]we briefly 
introduce the way to construct noisy weak lensing conver- 
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gence fields using N-body simulation. Measures of non- 
Gaussianity, i.e., the measures of the Gaussianization 
performance in the noisy case are presented in Sec lIIII To 
reduce the degradation of the Gaussianization method 
caused by shape measurement noise, we propose the 
Wiener filter as a remedy and redo the analysis in Sec lIVI 
In Sec|V] we summarize the results and outline the key 
issues for further investigation. 



II. SIMULATED LENSING CONVERGENCE 
MAPS WITH SHAPE MEASUREMENT NOISE 

We refer the readers to paper I (Yu et aL[l]) and the 
references therein for weak lensing basics. In the first 
step, we construct noise-free lensing convergence maps, 
as in paper I. Our N-body simulation was run using the 
Gadget-2 code [9|, adopting the standard ACDM cos- 
mology, with — 0.266, fl\ = 1 — n,n, as = 0.801, 
h = 0.71 and Ug = 0.963. The output redshifts are spec- 
ified such that any two adjacent outputs are separated 
by the box size L = 300/t~^Mpc in comoving distance. 
More details can be found in [l3|. We stack eight ran- 
domly shifted and rotated snapshots to comoving dis- 
tance 2400/i~^Mpc to make k maps. These lensing maps 
correspond to source galaxy redshift Zs ~ 1.02, with map 
size being 7.64° x 7.64°. When analyzing these maps, we 
use 512^ uniform grids, corresponding to a pixel size of 
0.9 . At a typical lens redshift zl = 0.5 of a typical source 
redshift Zs = 1, this corresponds to a size of 0.4/i~^Mpc, 
well into the nonlinear regime. This pixel size is suitable 
for investigating the nonlinear and non-Gaussian regime 
and for quantifying the performance of Gaussianization. 
We caution that the Gaussian transformation is nonlin- 
ear and its dependence on pixel size is nontrivial. We 
leave this issue for future work. 

In observation, the lensing convergence map con- 
structed from cosmic shear suffers galaxy shape measure- 
ment noise, with an rms of (Ttv — 0.2/^/N, where N is the 
mean number of galaxies per pixel. In the limit N — >■ oo, 
the shape noise distribution becomes Gaussian according 
to the central limit theorem. But for typical lensing sur- 
veys such as DES and LSST, N 10-40 and the shape 
noise may not fully reach a Gaussian distribution. Nev- 
ertheless, full investigation of the form that shape noise 
takes is well outside the scope of this paper. To proceed, 
we approximate the shape noise as Gaussian white noise 
for testing the Gaussianization method. We add two lev- 
els of Gaussian white noise into the convergence maps 
constructed from simulation, which correspond to source 
galaxy number densities of 10/arcniin^ and 40/arcniin^. 
Hereafter we denote the two cases as case A and case B. 
The first corresponds to DES-like survevs ITtII and the 
second corresponds to LSST- like surveys [l^- In real- 
ity, these surveys not only have different source galaxy 
number densities, but also different source galaxy red- 
shift distributions and hence they have different lensing 
signals. Since our purpose is to investigate the depen- 



dence of the Gaussianization on measurement noise, we 
neglect the difference in source redshift distribution and 
fix the lensing signal, but vary the measurement noise. 
For each noise level, we produce 20 Gaussian white noise 
maps and add them into convergence maps pixel by pixel. 
Throughout the paper, we denote the true lensing signal 
as kt, the noise as kn and the sum as k = kt + kn. 

For the given pixel size, the corresponding noise rms 
is ^ 0.07 for case A and ~ 0.035 for case B. The lens- 
ing signal rms is ~ 0.02. Hence, for both cases the noise 
dominates over the signal and significantly changes the 
one-point PDF, and pushes it towards a Gaussian form. 
The numerical results of one convergence map with two 
levels of noise added are presented in the left column 
of FiglH Indeed, for case A, there is almost no visible 
deviation in the one-point PDF from the Gaussian dis- 
tribution. For case B, only small deviations can be seen. 



III. GAUSSIANIZATION PERFORMANCE ON 
NOISY MAPS 

Here we briefly introduce the Gaussianization we use. 
Details are given in paper I. The Gaussian transformation 
K — y is deflned as follows: 



with normalization at the zero point, 
dy 



(1) 
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Here by definition Paiv) = exp(— j/^/2cr^)/v27ra-y is the 
Gaussian PDF and ay is the rms dispersion of y. We ob- 
tain this transformation numerically for each pixelized 
simulated map. The normalization coefficient of the 
K ^ y relation does not infiuence our measures of non- 
Gaussianity defined below, nor does it alter the perfor- 
mance of Gaussianization. The particular normalization 
both in paper I and here (Eq. [2]) is chosen such that 
(i) in the limit that the k field becomes Gaussian, the 
Gaussianization transformation becomes an unity trans- 
formation, and (ii) it shares the same unity slope as the 
logarithmic transformation at k = 0. 

By definition, the one-point PDF of the y field of the 
given pixel size is Gaussian. However, this does not nec- 
essarily mean that the resulting field y is a multivariate 
Gaussian random field. For example, if we smooth the 
field, the one-point PDF of the smoothed field ys can be 
non-Gaussian. Following paper I, we will evaluate vari- 
ous non-Gaussianity statistics, such as the cumulants as 
a function of smoothing scale and the bispectrum. The 
cumulants Kn {n = 1, 2, • • •) are defined by 
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FIG. 1: The left column shows the PDFs of one convergence map we constructed with two levels of noise added. The right 
column shows the PDFs of the convergence map with noise reduction by the Wiener filter. The PDFs of the field after the 
Gaussianization are also plotted (dotted line), which is just the standard Gaussian distribution by definition. From the PDFs we 
can see that the Gaussian noises of both levels dominate the distribution. After the Wiener filter is applied, the non-Gaussianity 
is recovered to some extent. 
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The reduced bispectrum is defined by 
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A. The K-y relation of noisy maps 



The Gaussian transformations we obtain are presented 
in the left panel of FiglS] For case A, the noise dominates 
over the signal {\kn\ > I'^tI) for almost all of the pixels. 
Since the noise PDF is Gaussian to a good approxima- 
tion, we have y = k. This is what we see from Fig. [51 
On the other hand, visible deviation from y = k can be 
found at the high k end, where the lensing signal can 
dominate over the noise due to its skewed PDF. 



with i'l + £2 + ^3 = 0. We express q in the coordinate 
a), where I = li, u = I2/I1 and a — arccos(Zi • 

l2)/hl2. 

Again, we caution the readers that the definitions of 
the cumulants and the reduced bispectrum are different 
from the widely used definitions in large scale structure 
study. They are designed to be invariant under the trans- 
formation y — Ay, where A is an arbitrary constant. 
This property is desirable when comparing these non- 
Gaussianity measures before and after the Gaussianiza- 
tion. 



For case B, noise is reduced by a factor of 2, so the 
non-Gaussianity of the PDF is less affected. As a conse- 
quence, we find larger deviation from y = k. 



The Gaussian transformations for 20 realizations con- 
verge with each other very well. The convergence is much 
better than what we found in the previous work, in which 
the lensing maps are free of noise. This better conver- 
gence is also a manifestation of overwhelming shape mea- 
surement noise. 
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FIG. 2: The left panel shows the Gaussian transformations from the noisy k to Gaussian y field, for 20 realizations. The right 
panel shows the result after the noise is reduced by the Wiener filter. The heavy noise makes the Gaussian transformations of 
trivial-mapping {y = n). Only a few pixels survive in case A. For the noise reduced case, we can see the effect of Wiener filter 
from the shape of the Gaussian transformation. The 20 Gaussian transformations converge with each other for noisy maps, 
and are in reasonable agreement with each other for the noise-reduced case. The small divergence in the latter case arises from 
the different normalization factors determined at the zero point. Nevertheless, the normalization factor will not influence the 
non-Gaussianity measures. 



B. The cumulants of smoothed noisy maps 

The cumulants of the smoothed filed Kn{ys) £^re non- 
trivial checks of non-Gaussianity, despite the fact that the 
cumulants of the unsmoothed field Kn{y) are Gaussian 
by definition. Smoothing introduces pixel-pixel correla- 
tion and gives rise to non-Gaussian ys cumulants. The 
results for the cumulants of the smoothed field up to the 
6th-order are presented in the left column of Figsl2]and|3] 
for case A, and in Figsl5]and[n]for case B. As in paper I, 
all results are averaged over 20 realizations and the error 
bars are the rms dispersions among these realizations. 
Note that these error bars on different smoothing scales 
are highly correlated with each other. To highlight the 
performance of the Gaussianization, Kn{ys)/ Kni^s) as 
functions of smoothing scale are shown in Figl?! 

Contamination of shape measurement noise signifi- 
cantly changes the behavior of the cumulants in two ways: 
(i) Since the noise field is Gaussian and the measure- 
ment noise often dominates the lensing signal, the cu- 
mulants Kn{Ks) are much lower than the corresponding 
ones in the noise-free case (paper I), (ii) Smoothing sup- 
presses the noise more efficiently than the signal, since 
the power of noise is more concentrated to small scales 



(e.g., zero-lag). Hence for sufficiently large smoothing 
scales, the signal can dominate over the noise and the 
non-Gaussianity of the signal can show up. This leads 
to larger Kn for larger smoothing scales for case A. On 
the other hand, smoothing decreases the non-Gaussianity 
of the signal. So for the noise-free case in paper I, Kn 
decreases with the smoothing scale. Case B has a noise 
level between case A and the noise-free case. The two ef- 
fects of smoothing compete and cause ^^4.5,6 to increase 
with the smoothing scale first and then decrease with it. 



The heavy shape measurement noise contamination 
significantly degrades the performance of the Gaussian- 
ization. For a shallower lensing survey like DES (case 
A), the Gaussianization has very limited capability to 
reduce non-Gaussianity (FigslSl |4]and[7l). As we expect, 
the Gaussianization works better for the case of lower 
noise. But even for deep surveys like LSST (case B), the 
Gaussianization can only moderately suppress the non- 
Gaussianity. It suppresses the skewness by a factor less 
than 2. Although the situation is better for the higher- 
order cumulants, the suppression factor is still less than 

10 for if4,5,6- 
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FIG. 3: The ensemble averaged skewness (top) and kurtosis (bottom) of the smoothed noisy (left) and the Wiener filtered (right) 
convergence maps are plotted in a solid line as a function of smoothing scale for case A. The results of the Gaussianized fields 
are plotted in a dotted line with a little horizontal shift for clarity. The error bars are the rms dispersions among 20 realizations. 
The Gaussianization cannot reduce the skewness and kurtosis of these noisy maps to an acceptable level. However, we find 
that the Wiener filter is efficient in reducing the shape measurement noise. For the Wiener filtered maps, the Gaussianization 
can suppress the skewness and kurtosis towards zero for all the smoothing scales we considered. The statistical errors are 
quite correlated, since K„ at different smoothing scales share much the same cosmic volume, and hence much the same cosmic 
variance and measurement noise. 



C. The bispectrum of noisy maps 

Another measure of non-Gaussianity is the reduced bis- 
peetrum q{li, h, h) = u, a), whieh vanishes for Gaus- 
sian fields. The reduced bispectrum for various configu- 
rations are presented in Fig[8] The Gaussianization al- 
most entirely fails to reduce the bispectrum for case A, 
in which the shape measurement noise dominates. So 
we only show the results for case B. To better demon- 
strate the impact of shape noise, we choose the same 9 
configurations of I G [200, 6000) as in paper I. For most 
configurations of a), q can only be suppressed by 
the Gaussianization by a small factor. For many config- 
urations, the qs after the Gaussianization deviate signif- 
icantly from zero, and stay positive. This result differs 
from the noise-free case in paper I, in which the g's after 
the Gaussianization are consistent with and scatter about 
zero. Since the bispectrum is an unambiguous measure 
of non-Gaussianity, these nonzero results show the fail- 
ure of the Gaussianization. This poor performance of the 
Gaussianization is consistent with the results in the pre- 



vious section on the skewness, which is the bispectrum 
weighted over all configurations. 



IV. GAUSSIANIZATION PERFORMANCE ON 
WIENER FILTERED MAPS 

The way that Gaussianization works is that it rela- 
tively down-weights large k pixels (regions) which are 
mostly responsible for the non-Gaussianity. Hence to im- 
prove the performance of Gaussianization for noisy lens- 
ing maps, we shall apply a filter to down- weight regions 
of low S/N, so regions of high S/N can be highlighted. In 
this way the power of Gaussianization can still persist. 
For this purpose, we will adopt the Wiener filter. The 
Wiener filter is a filter in Fourier space. But statistically 
speaking, it has the desirable property of highlighting 
high S/N regions. It is also widely used in astronomy. 
In particular, it is shown by [ll| that the Wiener filter 
can improve the reconstruction of one-point (true) con- 
vergence PDFs from noisy lensing maps. 
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FIG. 4: The ensemble averaged 5th- (top) and 6th- (bottom) order cumulants of the smoothed noisy (left) and noise-reduced 
(right) convergence maps are plotted in a solid line as a function of smoothing scale for case A. The results of the Gaussianized 
fields are plotted in a dotted line with a little horizontal shift for clarity. The correlated error bars are the rms dispersions 
among the 20 realizations. The Gaussianization cannot reduce the 5th- and 6th-order cumulants of these noisy maps to an 
acceptable level. When the noise is reduced by the Wiener filter, the Gaussianization can suppress the 5th- and 6th-order 
cumulants to zero for all the smoothing scales we considered. 



A. Wiener filter 

In signal processing, the Wiener filter is widely used 
to reduce the amount of noise presenting in a signal by 
comparison with an estimation of the desired noiseless 
signal. For simplicity, we briefly review the Wiener filter 
in ID space. Assume that the corrupted signal is s{x) — 
s{x) + n(x), in which s{x) is the original signal and n{x) 
is the noise. We can convolve s{x) with a filter h{x) to 
reconstruct the signal s{x). The reconstruction error has 
a dispersion 

E{H{k)) = y" dfc \H{k)[S{k) + N{k)] - S{k)\^ , (5) 

in which H{k), S{k), N{k) are the Fourier transforms of 
h{x), s{x) and n{x), respectively. The Wiener filter is 
the one to minimize E, 

H(k) 1^(^)1' (6) 

''^''>- \Sik)\^ + \NikW ■ ^ ' 

At scales where the signal dominates over the noise, the 
filter is close to 1, while at scales where the noise dom- 



inates over the signal, the filter is close to 0. This is a 
desirable property for our purpose. 

In our application, the |5(fc)p term in the Wiener filter 
is the power spectrum of the k field P{£) = {\k{£)\'^). The 
|A^(fc)P term is the Gaussian white shape noise power 
spectrum, which does not vary with scale. 

We will basically redo the analysis from Sec lIIII for 
noise reduced maps by the Wiener filter Eq.®. Since in 
observation the noise power spectrum will not be known 
exactly, we also test the Wiener filtering in which the 
noise power spectrum is overestimated or underestimated 
by a factor of 2 for case B. We may expect that the 
noise power spectrum estimated in a real survey would 
fall somewhere between the two extreme cases. 

The PDFs of the Wiener filtered convergence maps 
are presented in the right column of FiglTJ Compared 
to the unfiltered maps in the left column, the Wiener 
filter indeed recovers the lensing signal to some extent, 
and as a consequence the PDFs are now visibly non- 
Gaussian, even for case A, which has larger measure- 
ment noise. Nevertheless, the Wiener filter significantly 
improves the PDF recovery. Later we will show the gain 
of non-Gaussianity of the maps. But compared to the 
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FIG. 5: The ensemble averaged skewness (top) and kurtosis (bottom) of the smoothed noisy (left) and noise-reduced (right) 
convergence maps are plotted in a solid line as a function of smoothing scale for case B. The results of the Gaussianized fields 
are plotted in a dotted line. The correlated error bars are the rms dispersions among the 20 realizations. The Gaussianization 
cannot reduce the skewness and kurtosis of these noisy maps to an acceptable level. When the noise is reduced by the Wiener 
filter, the Gaussianization can suppress the kurtosis towards zero for all the smoothing scales we considered. But for the 
skewness, we still see residual non-Gaussianity. The results for an inaccurate noise spectrum estimation used in the Wiener 
filtering are presented as dashed and dot-dashed lines. A factor of 2 overestimation or underestimation does not change the 
Gaussianization performance much. 



strong non-Gaussianity existing in the noise-free case, it 
is obvious that the Wiener filter can not recover the field 
too much. For case A, the PDF is still close to a Gaus- 
sian distribution. The recovery effect is more obvious for 
case B. 



B. The K-y relation of noise- reduced maps 

The Gaussian transformations for the noise-reduced 
convergence maps are presented in the right panel of 
FiglH Compared to the left panel, the shape of the 
transformation k — > y apparently changes. The turning 
points between the straight line y = k (which represents 
the dominance of Gaussian noise) and the high k end 
(where the Gaussianization takes effect) shift towards a 
low K value. Though the reduced noise still makes the 
Gaussian transformation a trivial mapping for the low 
K value pixels, pixels with higher k persist through the 
noise thanks to the filtering, and thus the Gaussianiza- 
tion method will have effects on these high k value pixels. 



The 20 Gaussian transformations after noise reduction 
are in reasonable agreement with each other for both case 
A and B. The small divergence in this noise- reduced case 
rises from the different normalization factors determined 
at the zero point. Nevertheless, the normalization factor 
will not influence the non-Gaussianity measures. 



C. The cumulants of noise-reduced maps 

The cumulants 1^3.4,5.6 for the noise-reduced case are 
presented in the right column of Figs El and |3] for case 
A and Figs [5] and H] for case B. The solid lines are for 
the filtered convergence maps prior to the Gaussian- 
ization KninY^), while the dotted lines are for after- 
wards Kniug^), ill which superscript WF means that 
the noise is reduced by the Wiener filter. All the re- 
sults are averaged over 20 realizations and the error bars 
are the rms dispersions. The ratios of the cumulants 
Knius'^) / Kn{Kg^) presented in the right column of 
FigiZl 
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FIG. 6: The ensemble averaged 5th- (top) and 6th- (bottom) order cumulants of the smoothed noisy (left) and noise-reduced 
(right) convergence maps are plotted in a solid line, as a function of smoothing scale for case B. The results of the Gaussianized 
fields are plotted in a dotted line. The correlated error bars are the rms dispersions among the 20 realizations. The Gaussian- 
ization cannot reduce the 5th- and 6th-order cumulants of these noisy maps to an acceptable level. When the noise is reduced 
by the Wiener filter, the Gaussianization can suppress the 5th- and 6th-order cumulants to zero for all the smoothing scales we 
considered. The results for an inaccurate noise spectrum estimation used in the Wiener filtering are presented as dashed and 
dot-dashed lines. A factor of 2 overestimation or underestimation does not change the Gaussianization performance much. 



The maps with noise reduction have more skewness 
than the case without noise reduction (left column) , and 
have less skewness than the noise-free case in the former 
work, which is just as we would expect from the PDF re- 
sult. On large smoothing scales, the recovery of the non- 
Gaussianity due to the smoothing out of the noise and the 
reduction of the non-Gaussianity due to the smoothing 
of the signal, seem to just cancel each other for case A, 
leading to an approximately constant skewness. Similar 
results also occur for K^^^ q. In case B, for -fC4,5.6 on large 
smoothing scales, the reduction effect dominates over the 
recovery effect. We can see the cumulants stay constant 
first on small smoothing scales due to the equilibrium of 
the two opposite effects, and decrease on large smoothing 
scales due to the overwhelming reduction effect. 

For both case A and B, for all order cumulants we con- 
sidered, the Gaussianization works much better than the 
case without noise reduction. But we still see nonzero 
skewness on the largest smoothing scale. The kurtosis 
and the 5th- and 6th-order cumulants of the Gaussian- 
ized maps are well consistent with zero within errors. The 
Gaussianization also reduces the variance among differ- 



ent realizations. 

We can quantify the performance of the Gaussianiza- 
tion in the right column of Figd For both case A and 
B, the Gaussianization can suppress the skewness by a 
factor of 5 on the largest smoothing scale, and better 
(more than 10) in higher-order cumulants. The Gaus- 
sianization has much better performance on the Wiener 
filtered case (as we expected), since the Wiener filter re- 
duces the noise and recovers the signal to some extent, 
and the Gaussianization has an effect on more pixels that 
are dominated by the non-Gaussian signal. 



D. The bispectrum of noise-reduced maps 

The reduced bispectrum of noise-reduced maps for var- 
ious configurations are presented in Fig 121 where the solid 
lines are for the noise-reduced fields q{K^^) and dotted 
lines are for the corresponding q{y^^) fields. The results 
are averaged over the 20 realizations and the error bars 
are the rms dispersions within them. We show the bis- 
pectrum on the same 9 configurations as the g's in the 
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FIG. 7: The ratios of the nth-order cumulants of the smoothed noisy (left) and noise-reduced (right) convergence maps from 
prior to the Gaussianization to afterwards, as functions of smoothing scale for case A (top) and B (bottom). The Gaussianization 
can only suppress by a factor of 2 for the skewness of a noisy map, and a little more or less for higher order cumulants. When 
the noise is reduced by the Wiener filter, the Gaussianization can suppress by a factor of 5 for the skewness of a noisy map, 
and more than 10 for higher-order cumulants. 



last section and in the former work. For most configu- 
rations of (/, u, a), prior to the Gaussianization 
have larger values than the results in the previous section, 
due to the recovery of non-Gaussianity by the Wiener fil- 
ter. We can see that after the noise reduction, for many 
configurations g's can be significantly suppressed by the 
Gaussianization, and the g(t/^^)'s are consistent with 
zero within error bars. These results show the effective- 
ness of the Gaussianization after the noise reduction by 
the Wiener filter. This conclusion is consistent with the 
findings in the previous section on the skewness. 

Since q describes the non-Gaussianity of an individ- 
ual configuration, to some extent it is a more sensitive 
measure of the performance of the Gaussianization than 
skewness. We find robust evidence for residual non- 
Gaussianity in some configurations of the bispectrum 
(e.g., some configurations with large a in Fig. |9]). This 
means that the non-Gaussianity and pixel correlation 
is also not completely removed by the Gaussianization, 
which is the same situation as in the former work. A di- 
rect implication is that the Wiener filter can not fully re- 
cover the copula of the convergence field. Another source 
is the small deviation from the Gaussian copula of the 
noise- free convergence field, as we found in paper I 



E. Inaccurate noise estimation in the Wiener filter 

We test whether a precisely known noise power spec- 
trum is crucial in Wiener filtering. We mimic uncertain- 
ties in noise power spectrum determination by artificially 
multiplying the correct noise power spectrum by a factor 
of 2 or 1/2. The results are presented in the dashed line 
and dot-dashed line in FigsEl El and O for case B. The 
Gaussianization is indeed affected by uncertainties in the 
noise power spectrum determination. However, even for 
a rather large overestimation or underestimation by a 
factor of 2, the Gaussianization remains efficient. This 
behavior is more or less what we expect. The Wiener 
filter works for the Gaussianization since it approaches 
one where signal overwhelms noise, and approaches zero 
where noise overwhelms signal. A wrong noise power 
spectrum does not alter this crucial property; it only af- 
fects the transition region where signal is comparable to 
noise. As long as the noise estimation is reasonably good, 
the Wiener filter would be robust in improving the Gaus- 
sianization performance. 

Furthermore, we notice an interesting fact: an inaccu- 
rate estimation of the noise spectrum used in the Wiener 
filter may improve the performance of Gaussianization. 
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FIG. 8: The ensemble averaged reduced bispectrum q{£, u, a) of the noisy convergence maps for case B on the same 9 configura- 
tions as the previous work. The qs prior to the Gaussianization are presented in a solid line, while the afterwards is in a dotted 
line with a horizontal shift for clarity. The error bars are the rms dispersions over 20 respective maps. The Gaussianization 
can only suppress the bispectrum by a small factor, and a large amount of non-Gaussianity still remains in bispectrum of the 
new field. 



This implies that the Wiener filter may not be the op- 
timal filter to improve the Gaussianization. This issue 
should be investigated further. 



V. DISCUSSIONS AND CONCLUSIONS 

We have demonstrated that the Gaussianization 
method we proposed in the former work fails when Gaus- 
sian shape noise is added into the convergence maps. 
We quantify the performance of the Gaussianization by 
various measures of the non-Gaussianity, such as skew- 
ness, kurtosis, the 5th- and 6th-order cumulants of the 
smoothed fields, and the bispectrum. We found that the 
Gaussianization can only suppress the skewness by a fac- 
tor less than 2, and more or less for higher-order cumu- 
lants for different noise levels. This implies that we can 
not treat the resulting y field as Gaussian, since higher- 
order statistics are still large due to the noise corruption. 

We propose the Wiener filter to improve the perfor- 
mance of the Gaussianization for the real noise existing 
case. We find that the Gaussianization works well on 



the Wiener filtered convergence maps. We also find that 
the performance is rather robust against uncertainties in 
determining the noise power spectrum in the Wiener fil- 
tering. So an accurate determination of the noise power 
spectrum is not crucial in applying the Wiener filter and 
the Gaussianization to observational data. This insensi- 
tivity improves the applicability of the Gaussianization 
to real data. The Gaussianization procedure can trans- 
port cosmological information in higher-order statistics 
of noisy lensing fields into the power spectrum of y. For 
this reason, analyzing weak lensing statistics can be sig- 
nificantly simplified. Nevertheless, measurement noise 
causes irreducible corruption in the Gaussianization per- 
formance, compared to the noise-free case. 

Due to the lack of fully independent realizations, finite 
box size and limited number of particles of the simu- 
lation, our findings are inevitably affected by numerical 
issues such as cosmic variance, rnass and force resolution, 
shot noise and aliasing effect [i2[ . We do not expect that 
any of these factors will alter the major conclusion on the 
failure of Gaussian transformation due to the heavy noise 
and the success after the noise reduction by the Wiener 
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FIG. 9: The ensemble averaged reduced bispectrum q{£, u, a) of the noise-reduced convergence maps for case B, on the same 
9 configurations. The qs prior to the Gaussianization are presented in a solid line and those afterwards in a dotted line with a 
horizontal shift for clarity. The error bars are the rms dispersions over 20 respective maps. The recovery effect of the Wiener 
filter results in a larger value in the bispectrum. The Gaussianization after noise reduction can suppress the bispectrum more, 
and on many configurations the bispectrum afterwards is close to zero within error bars. But there still exists some amount 
of non-Gaussianity on small scale configurations. The results for an inaccurate noise power spectrum estimation used in the 
Wiener filtering are presented in dashed and dot-dashed lines, for cases of a factor of 2 overestimation or underestimation 
in the noise power spectrum. The results show that, even for these extreme cases, the Gaussianization performance remains 
impressive. 



filter. These numerical issues could in particular have 
a larger impact on the residual non-Gaussianity, which 
itself is weak and relatively hard to measure accurately. 
However, it is very unlikely that the detected residual 
non-Gaussianity can be canceled exactly. In this sense, 
the detection of the residual non-Gaussianity is robust, 
although we need many more simulations to measure the 
amplitude to high precision. 

A number of works have addressed extra complexities 
in Gaussianization. Neyrinck et al.^ checked the effect 
of the logarithmic transformation and Gaussian transfor- 
mation on the galaxy density field. They found that — 
although due to the presence of discreteness noise the 
galaxy density field cannot be fully Gaussianized — both 
transformations still dramatically reduce nonlinearities 
in the power spectra of cosmological matter and galaxy 
density fields. But the transformations do increase the 
effective shot noise. Neyrinck 4] studied the sensitiv- 
ity of five cosmological parameters constrained to the 



Gaussianized power spectra by logarithmic transforma- 
tion and the Gaussian transformation. He found the 
power spectrum of the log-density provides the tightest 
cosmological parameter error bars in all five parameters 
tested. From our studies, we found that for both the 
noise-free case and the noisy case the Gaussian transfor- 
mation apparently deviates from the logarithmic trans- 
formation. Although the logarithmic transformation will 
not produce a perfect Gaussian random field, there is no 
free parameter in the transformation. Thus the informa- 
tion will not be lost into the freedom of the transforma- 
tion form, and theoretical predictions can be made, such 
as [la] , 

Joachimi et al.^ employed an extended form of the 
logarithmic transformation, using Box-Cox transforma- 
tions k{X, a) — [{k + a)'^ — 1]/A with two free parameters 
(A 0, a) to Gaussianize ideal weak lensing convergence 
field. (The Box-Cox transformation reduces to the loga- 
rithmic transformation when A = 0.) The optimized Box- 
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Cox transformation is determined by fitting it to a Gaus- 
sian transformation through the free parameters. They 
developed analytical models for the transformed power 
spectrum, including effects of noise and smoothing, and 
found that the optimized Box-Cox transformation per- 
forms better than an offset logarithmic transformation 
in Gaussianizing the convergence; but, none of them are 
capable of eliminating correlations of the power spectra 
between different angular frequencies. When they add a 
realistic level of shape noise, all the transformations per- 
form poorly with little decorrelation of angular frequen- 
cies, and the arc-tangent logarithmic transformation — 
which approximates a straight line near the zero point — 
was proposed to deal with the shape noise situation. 

Carron [l3| shows a peculiarity of fields with heavy 
tails: the hierarchical statistics (N-point functions) do 
not entirely specify the distribution. Thus the cosmolog- 
ical information is not only pushed towards higher-order 
statistics due to the nonlinear evolution, but can also be- 
come inaccessible through the extraction of the full series 
of moments of the field. His more recent work [l^ shows 
that the entire hierarchy of moments quickly ceases to 
provide a complete description of the convergence 1-point 
PDF when leaving the linear regime. However, a simple 
logarithmic mapping makes the moment hierarchy well- 
suited again for parameter extraction. Carron fl^ shows 
that the power spectrum of the log-density field carries 
more information than the power spectrum of the field 
when entering the nonlinear regime. These works indi- 
cate that the logarithmic transformation and the Gaus- 
sianization method may work beyond simply propagating 
cosmological information in the moment hierarchy, and 
support the application of these methods in data analysis 
and cosmology. 

Nevertheless, there are many key issues to be explored 
in future studies. An incomplete list includes the follow- 
ing. 

• The applicability to real data. Although in this 
work we consider a more realistic case than in the 
former work, it still differs from real observations. 
The application is complicated by various measure- 
ment errors in real data. As long as the pixel size is 
sufficiently large, the central limit theorem drives 
its distribution to be Gaussian. Since the Gaus- 



sianization that we have proposed is nonlinear, it 
will render this Gaussian noise into a non-Gaussian 
one. The Wiener filter could make the situation 
more complicated. Even worse, the same nonlinear 
transformation and the Wiener filter could mix the 
lensing signal and measurement noise. 

• The residual non-Gaussianity. The proposed Gaus- 
sianization works expectedly worse in the noisy case 
and the Wiener filter deals with the noise surpris- 
ingly well, but it does not perfectly produce a Gaus- 
sian random field. The residual non-Gaussianity 
we have detected is weak and is unlikely to carry 
a significant amount of cosmological information. 
Nevertheless, it may still be worth investigating 
the cosmological information carried by the resid- 
ual non-Gaussianities. 

• Other methods dealing with noise corruption. The 
Gaussianization we proposed fails due to the noise 
corruption, and succeeds in the noise-free case. The 
results of the Wiener filter set in somewhere in be- 
tween. There may exist some other appropriate 
tool dealing with the noise corruption which mean- 
while keeps the effectiveness of the Gaussianization. 
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